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The flow of a thermoelectric current through a semiconductor of submicron dimensions is an- 
alyzed. The rate of surface relaxation of the energy is assumed to be much higher than the rate 
of electron-electron collisions. Under these conditions, it is incorrect to describe the electron gas 
by means of a Maxwellian distribution and thus to describe the thermoelectric effects in terms of 
an electron temperature and a chemical potential. A theory is derived for these effects. This the- 
ory does not include the latter parameters and is based on a non-Maxwellian distribution which is 
spatially nonuniform in terms of energy. 



In the conventional theories on thermoelectric effects, the temperature and the chemical potential are fundamental 
parameters used by the theoretician. As was shown in Ref. |l|, on the other hand, when a submicron-size semiconductor 
layer makes contact with a heater and a refrigerator, the symmetric part of the distribution function of the current 
carriers is essentially non-Maxwellian (non- Fermi) , so such concepts as temperatures and chemical potentials are not 
at all valido There is accordingly a need to derive a theory of thermoelectric effects which is not based on these 
concepts. In the present paper we attempt to derive such a theory. 

For definiteness we consider a model of a closed thermoelectric circuit consisting of a semiconductor layer of thickness 
2a, an external metallic region of lw&'th L, and transition regions (contacts) with thicknesses 26 — > between the 
semiconducting and metallic mediaO'El We choose the thickness 2a to be so small that the inequality 2a <C l e holds 
(we are thus dealing with a submicron-size layer), where l £ , is the diffusion length for the relaxation of the energy of 
electrons among scattering centers. For quasielastic scattering events, this length is on the order of 1-10 fim. The plate 
surfaces x = =F(a + 5) are in contact with heat reservoirs at temperatures T\ and T 2 (Tj > T2). As a result, there are 
spatial variations in the distribution function, and a voltage is generated in the circuit. The impact of these variations 
increases with increasing frequency of the surface energy relaxation, v s . Since the latter is inversely proportional 
to the layer thicknesses^ the rate of electron-electron collisions in a submicron layer at a fixed equilibrium density 
may be much lower than v s . The electron gas would therefore not be correctly described by means of a Maxwellian 
distribution function with an electron temperature T e . In this case the symmetric part of the distribution function 
should be determined from a Boltzmann kinetic equation. 

As the mechanism for the inelastic interaction of the nondegenerate electron gas with the heat reservoirs we adopt a 
mechanism involving the flow of a "partial" current, by which we mean a current of a group of electrons with a given 
energy, through surface layers. Obviously, energy may be dissipated both at centers within the surface layers and in 
the course of a direct interaction with the heat reservoirs. Since our purpose here is not to derive exact quantitative 
results, but instead to establish fundamental positions ofithis theory, the energy dissipation within the surface layers 
is not of fundamental importance and it can be ignored.UI 

The symmetric part of the electron distribution function calculated uiwier these assumptions, in the approximation 
linear in the parameter AT/T < 1 [AT = T 1 -T 2 ,T= (7\ - T 2 )/2], isllJ 
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jo = j / {AT/T); j is the thermoelectric current density; a — J °° decr(e) is the bulk electrical conduct ivity: e, m, and 
£ are the charge, effective mass, and energy of an electron; /1 is the chemical potential at the temperatureQ T; fj,(e) is 
the density of electron states; t(s) is the momentum relaxation time; and 



Cx = — f°de£(e)#(e) 



and C2 are constants to be determined. The subscript s refers to surface layers. 

It follows from these definitions that the functions <j(e) and <j s (e) constitute partial bulk and surface conductivities. 
We have £(e) = (T/e)a(e). For a surface layer, we cannot, in general, write a similar relation for £ s (e) and cr s (e) 
at the outset, since .the kinetic coefficients of the boundary layer change with decreasing thickness 2<5, and they may 
acquire singularities^ in the limit 28 — > 0. It is not possible to follow these changes in the absence of a microscopic 
theory for surface relaxation processes. We will accordingly assume for simplicity £ s (e) = (T/e)/3 s a s (s), where (3 S is 
a phenomenological surface parameter with the dimensionality of a reciprocal length. We do not believe that this 
simplification is of fundamental importance. The important point is that there is no relationship of any sort between 
the functions a(e) and cr s (e). 

Under this assumption we can rewrite the functions in \&(e) as follows: 
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The first term in Eq. (|^) corresponds to the effect of the heat reservoirs on the distribution function in the case in 
which these reservoirs are related to the electron gas of the submicron layer by the part of the partial currents caused 
by the internal thermoelectric field and by thermal diffusion [in the temperature approximation, ^(e) = M(s) — 
3/4 — e/2T; Ref. [5J]. The second term reflects the mechanism for the formation of a non-Maxwellian distribution 
function involving the flow of a macroscopic current across the interface between two media differing in the energy 
dependence of the relaxation time.EJ The last two terms stem from the drift components of the partial current in the 
semiconductor and in the transition layer. 

To find the constants C\ and C2, we use the definition of the coefficient C\ and the normalization condition on the 
function *(e) (Ref. @): 

/.DO 

de^{e)e- £/T - 0. (5) 







As a result, we find the system of algebraic equations 

71 Ci + 72C2 = -73 - 74 jo, 



7 1 Ci-7 2 C 2 = 73 + 74.70. 
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are constant, where 77(e) = ^fee~ £ l T . 

Determining the coefficients C\ and C2 from system of Eqs. (g), and substituting the results into Eq. (||) and ([!]), 
we find an expression for the symmetric part of the electron distribution function in the submicron layer: 



f (s,x) = e^' T jl + X - [F(s) + G(e)j ] 



(8) 



where 

F(e) = A(e)M(e) + r x Q( e ) + T 3 R(e), (9) 
G(e) = P(e) + r 2 Q(£)+r 4 i?(£), (10) 
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The thermal emf generated in the closed circuit is found from the expression! 

E = jR=J>?-dx = Km([ - dx + [ - dx + [ ^ dx + [ - dx V (12) 
/ cr stoo yj_ a+s a J a _ s a J a+S a J_ a _ s a J 

Here R = R s + R c + Rl is the total ohmic resistance of the closed circuit, R s = 2a/ a, is the resistance of the surface 
layer, R c = \im5 to o[2S / J °° dea s (e)] is the resistance of the contacts, and Rl = L/a^ is the resistance of the external 
part of the circuit, whose conductivity a^. 

For simplicity we set the cross-sectional area of the contour equal to one everywhere. 

The first and third integrals in Eq. ( |l2| ) correspond to the components of the thermal emf due to the semiconductor 
and the external region. The second and fourth integrals determine the thermal emf in the contacts. 

The electric current density j can be expressed in terms of the symmetric part of the distribution function: 



dej(e,x), (13) 







where the partial current j(e, x) is 



3m \ dx de 
E = —dtp/dx is the electric field, and ip is the electric potential 



^) = -^^f^+ei#), (14) 
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If we ignore the contribution to the thermal emf in the metallic part of the circuit, the integration over the closed 
contour in Eq. ( |l2|) reduces to 



2AT 



/ de<r(e)F(e) ~ j — / dea(e)G(e) 
Jo ecr A) 

dea s (e)eS £ ^^ /T [foo(e,x = a) - f a (e,x = a) - foo(e,x = -a) + f (s,x = -a)] . (15) 



jR = - : 



Here foo(e,x = ±a) are Fermi distribution functions for the electrons in the metal, with temperatures T2 and T%, 
respectively. 

Setting foots, x = a) « foo(s,x = — o), and carrying out some simple calculations, we find an expression for the 
thermoelectric current j: 
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Since the function A(e) = a(e)/a — a s (e)/a s appears linearly in G(e) [see Eq. ( |To| ) and (|ll|)], the denominator of 
Eq. ( |l6| ) has no singularities. 

After we substitute Eq. ([jl]) into Eq. (f[2|). we can write the thermal emf in the customary form 



E = aAT 



(17) 



where the thermal-emf coefficient a is 
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If the electron gas in the semiconductor is Maxwellized, the expression for the thermal-emf coefficient become: 

(1 - (3)a s + /3a v 
a = ^ ^ K, 



si 



R + Rp 



(19) 



where the coefficient determines the extent to which the thermal contact is isothermal, a v and a s are the bulk and 
surface components of the thermal emf, respectively, and Rp is the Peltier resistance. 
It can be concluded from a comparison of Eq. (|l|) and (]l|) that the quantity 



R>=?L 
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ds 
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G(e) 



is the Peltier resistance for the conduction electrons described by distribution Eq. (|I|). If <j(e) = cr s (e), then the 
thermal current and thermal emf vanish according to Eqs. (|l^) and (|l"8|). The evident reason for this vanishing is 
that the bulk and surface components of the overall thermal emf cancel out under the assumption that the partial 
conductivities tr(e) and cr s (e) are proportional to, respectively, the parameters £(e) and £ s (f), which are responsible 
for the thermal coupling of the electrons of the submicron layer with the heat reservoirs.El In this case, the Peltier 
resistance clearly must vanish automatically. 

Let us consider two particular cases corresponding to extreme values of the function Z(e): 2(e) <C 1 and 2(e) 3> 1. 
It is easy to show, on the basis of Eqs. (B|) and (0), that we have F(e)toQ in the first case: this circumstance 
means that there is no thermal current in the circuit. In this case, distribution function Eq. (|l|) transforms into a 
Maxwellian equilibrium distribution function with an equilibrium temperature T. This situation could naturally be 
called "adiabatic." Rewriting the inequality 2(e) < las a(e)/a ^> /3 s a s (e), we see that, for given values of the surface 
parameters, the "adiabaticity" of the system becomes more pronounced with thinner samples, i.e., with increasing 
partial conductivity of the semiconductor layer per unit length. This situation is equivalent to an effective decrease 
in the surface conductivity and thus a weakening of the coupling of the carriers with the heat reservoirs. In the 
temperature approximation, this situation corresponds to the condition that the thermal conductivity per unit length 
is greater than the surface thermal conductivity.!] In this case, the temperature T e , becomes uniform along the length 
of the sample because of the large bulk thermal conductivity. 
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In the second limiting case [Z(e) ^> 1] we have F(s) = M(s), G(e) — > 0, and the thermal-emf coefhcient is 

±-J™ decree) 



T 



(20) 



where the number r determines the relaxation time of the momentum of the electrons in the semiconductor: t(e) 
Me/Ty. 

In this case the distribution function (TH) becomes a Maxwellian distribution corresponding to the temperature 



x AT 

TJx) = T . 

y ' 2a T 



(21) 



This result may look paradoxical, as contradicting our original assumption that it is incorrect to describe the 
electrons by means of a Maxwellian distribution function. In speaking of a "Maxwellization" of a gas of carriers one 
is usually assuming that a dominant role is played by an electron-electron collision integral; such an integral can 
be omitted from the problem at hand. Consequently, the conventional mechanism for the shaping of a Maxwellian 
distribution function is actually meaningless. In submicron layers, however, one can point out another mechanism for 
a Maxwellization of a gas: the intense "mixing" of the partial electron current densities as a result of their interaction 
with the heat reservoirs. The heat reservoirs, in fact, determine the temperature (|2l]) of the electron gas. Clearly, as 
the frequency v v , decreases, this mixing weakens, and the distribution function becomes non-Maxwellian. It is also 
clear that the boundary conditions should be quite different from adiabatic conditions in this case, since otherwise 
the thermoelectric effects would vanish. The boundary conditions describing the second limiting case can be found as 
isothermal boundary conditions which lead to an ideal coupling of the electrons with the heat reservoirs. 

In the temperature approximation, the thermal emf is usually found fromlij 



E=l(-t 
e \ T 



I I AT 



(22) 



It contains a thermal-diffusion component Et — (l/e)(r + \)AT and a contribution from the contact potential 
difference, Ec — (l/e)(3/2 — /i/T)AT. Under nonisothermal boundary conditions, a third component is added: the 
surface thermal emf in ( |l9| ) . This partitioning of the thermal emf into components can be carried out because physical 
quantities — the chemical potential and the temperature — are distinguished in the Maxwellian distribution function. 
The temperature plays a twofold role in thermoelectric phenomena. First, the temperature gradient determines the 
"external— force" acting on the carriers. Second, the same gradient specifies the spatial variation of the chemical 
potentialLlJ and thus the internal thermoelectric field. 

Comparing Eq. (^0|) and (]2^), we see that in submicron semiconductors under isothermal boundary conditions 
it is again possible to separate contributions from thermal diffusion atid contact components of the thermal emf. 
Furthermore, the contact emf can be determined from an expression which includes exclusively surface characteristics: 



En = - 



3 _ j-J™deeo s {e) 
2 T 



AT. 



(23) 



Under nonisothermal boundary conditions, it is no longer possible to partition the thermal emf into bulk, contact, 
and surface components, as can be seen from expression (|lq). All we can do is speak in terms of "bulk" and "surface" 
components of the thermal emf, without subdividing the latter into "contact" and "thermal-diffusion" parts. 

The two limiting cases discussed above can occur only if the inequalities written above for the function Z(e) hold 
for all or at least most of the electrons. If the energy dependence of the relaxation time in the interior is different 
from that in the surface layer, then (for example) adiabatic boundary conditions may hold for the electrons of one 
group, and isothermal boundary conditions for those of another. The behavior of the overall electron gas will then be 
determined by an effective resultant contribution of each of these groups. 
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